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1 Introduction 



The region of parameter space for which the electroweak baryogenesis scenario is viable in the con- 
text of the Minimal Supersymmetric Standard Model (MSSM) is currently being tested at existing 
accelerators. Electroweak baryogenesis requires a strong first-order electroweak phase transition. 
This ensures that sphaleron transitions are switched off in the broken phase, permitting the preser- 
vation of a baryon asymmetry [[[], 0] (for reviews, see @-f|). In the context of the Standard Model 
the only free parameter relevant to the phase transition is the Higgs mass. Non-perturbative results 
of the phase transition show that for the Standard Model there is no experimentally allowed value 
of the Higgs mass for which the phase transition is sufficiently first order JF], |J0. Non-perturbative 
studies of the phase transition are necessary at finite temperature due to the existence of infrared 
divergences in gauge theories. However, recent lattice simulations have shown that for the case 
of the MSSM the perturbative results are conservative in the constraints they impose on the pa- 
rameters 0. Perturbative analyses of the phase transition have pointed out different mechanisms 
that can enhance the strength of the phase transition in the MSSM with respect to the Standard 
Model case. The fact that there are additional scalar particles in the spectrum, which strongly 
couple to Standard Model fields, thus contributing significantly to the finite-temperature effective 
potential opens some room for electroweak baryogenesis. The main conclusion from the perturba- 
tive analysis in the 4D theory [p!0||-[l^| was that low values of the ratio of the vacuum expectation 
values of the two Higgs doublets tan/3 = ^, large values of the pseudoscalar mass tua and a very 
light right-handed stop were favoured. Given the current experimental limits on the masses in this 
model incorporating two-loop corrections QCD corrections is necessary in order to determine the 
allowed parameter space for which electroweak baryogenesis can take place fL6| |- |19| . In addition 
the presence of the trilinear terms in the potential is necessary in order to provide the required 
CP-violating sources for the production mechanism of electroweak baryogenesis. We recall that the 
ratio of the vacuum expectation value of the scalar field to the temperatyre, determines the rate 
of the sphaleron transitions in the broken phase and non-zero values for the trilinear couplings tend 
to reduce the strength of the phase transition in the Higgs doublet direction if the other parame- 



ters are kept fixed |Ij§, ED]. This is because the mixing in the squark sector reduces the light stop 



2 Non-perturbative analyses also show good agreement with perturbation theory for Higgs mass values mh ;$ 70GeV 
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contribution to the cubic term in the potential, it increases the critical temperature of the phase 
transition and radiative corrections make the Higgs heavier, thus weakening the strength of the 
phase transition. However, ref. showed that the upper bound on the Higgs mass was arising 
from a zero-temperature effect and not from the requirement of having out-of-equilibrium sphaleron 
transitions. A lighter stop can compensate for the three effects that weaken the phase transition. 

Alternatively the analysis of the phase transition can be performed in the context of a bosonic 
effective 3D theory obtained using dimensional reduction [0, fl2"l~l - Perturbative calculations 
within the effective theory reproduce in a relatively simple way the results obtained with the 4D 
perturbative analysis. The additional benefit of the construction of the purely bosonic effective 
theory is that it provides a link between the parameters employed in 3D lattice simulations with 
the underlying 4D parameters of the theory. The parameters in the 3D Lagrangian are obtained 
using dimensional reduction at high temperature by matching the static Green's functions in the two 
theories, to a given order in the perturbative expansion, by integrating out the non-zero Matsubara 
modes with masses of the order of irT, where T is the temperature. A further reduction can also 
be performed noting that some of the static modes in the theory have acquired thermal masses 
proportional to a gauge coupling multiplied by the temperature, ~ g w T, g s T. These so-called heavy 
particles can then be integrated out as well. References 0, 0, J26|, PTfl- IPOfl , give more details 
concerning the construction of effective theories for both the Standard Model and the MSSM. 

Interestingly the analysis of the phase transition with a light stop shows that for some region of 
parameter space a possible two-stage phase transition, first into a colour and and charge breaking 
(CCB) minimum and subsequently to the broken SU(2) minimum occurs. Initial lattice analysis 
suggest that the second transition is too strong and might not have taken place on cosmological 



time scales @. Furthermore, in a recent paper Cline et al. []3Ifl showed that although tunneling into 
the CCB minimum is viable the next step of getting out of this minimum appears not to be possible. 
Given this result, an accurate determination of the critical temperatures of the phase transitions 
can exclude an additional region of parameter space. Here we closely follow the procedure of refs. 
i n order to determine very precisely the critical temperature of the phase transitions. In 
particular, the 3D scalar masses require ultraviolet renormalization and a two-loop calculation (in 
4D) must be performed to fix the scales entering the mass parameters. This requires a full 4D 
effective-potential calculation incorporating mixing effects arising from the trilinear terms in the 
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MSSM potential. In this paper we extend the analysis of ref. [BU| to identify the range of parameter 



space for which a sufficiently strong first-order phase transition can occur and when a two-stage 
phase transition can exist. 

The perturbative component of the analysis relies on the validity of the high temperature ex- 
pansion. The value of the masses of the particles which are integrated out will define the regime of 
validity of the approach. It is interesting to consider the effects on the parameters of the effective 
theory as the mass of the third-generation left-handed squark doublet is varied. We present the 
formulae that allow us to continously pass from the low-mg ~ 300 GeV to large-mg ~ 1 TeV regime 
at one-loop with squark mixing and, for the latter regime, we compare our with results presented in 
the literature. In fact as we will discuss below for values the most relevant effects from the trilinear 
couplings are one- loop effects. For large values of tuq, the high temperature expansion is no longer 
valid; in this case we estimate the two-loop effects of the trilinear couplings. 

The paper is organized as follows: in section 2 we present the dimensional reduction to the 
effective bosonic theory at one-loop valid for values of itlq ~ 300 GeV. Section 2.2 presents a 
further one-loop reduction in the 3D theory, eliminating the heavy fields. In section 3 we present 
our results for the critical temperatures and the strength of the phase transition. The allowed region 
for electroweak baryogenesis to occur is also given here. Finally, in section4 we conclude. Appendix 
A presents the one-loop formulae needed when the trilinear terms are included and shows how the 
contributions to the 3D masses and couplings are modified for large values of third-generation left- 
handed squark doublet mass. In appendix B we give the expression for the two-loop unresummed 
effective potential in 4D, fully incorporating mixing effects, which is necessary for evaluating A# 3 , 
Ajj 3 . The contribution from the "heavy" particles that were integrated out at the second stageis 
given in this appendix. We discuss in appendix C the relevant zero-temeperature effects that must 
be included in our analysis. 

2 Dimensional Reduction 

We will now perform dimensional reduction by matching, as was previously done in refs. [26 



[f3(J for different models. Our initial 4D Lagrangian corresponds to the MSSM in the large-m^ 
limit. The particles that contribute to the thermal bath are the Standard Model particles plus 
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third-generation squarks: t L ,bL,tn,b R . There are two stages of reduction. The first one corresponds 
to the integration out of all non-zero Matsubara modes, that is with a thermal mass of the order of 
~ 7rT. We calculate all one-loop contributions to mass terms and coupling constants of the static 
fields to order g 4 , where g denotes a gauge or top Yukawa coupling. The second stage of reduction 
corresponds to the integration of heavy particles with masses of the order of g w T, g s T. 

2.1 First stage 

The potential in the 3D effective theory after integration over non-zero Matsubara modes is of the 
form 



V = m 2 H ,WH + X H3 (H^H) 2 + m 2 Ua U^U + \ U3 {U ] U) 2 + lz {H^ H){U ] U) 

+ rn 2 Qs Q^Q + m 2 D ,D^D + A®(H^H)(Q^Q) + A C 4 (H^Q)(Q^H) 

+ (AJ + hf^H^l 2 + (h? u + fflQi&QvUp 

+ tf?U a UZQ* jy Q jy + g?°D a D* Q* jfj Q ja 

+ gf 2 D D a D* a Q* n Q n + g v D U a U;D;D a 

+ g^UaU^D, + ArtQtQ) 2 + X Ds (^D) 2 + \ Q MQif 

+ df^QiaQjaQi-yQjl df^QiaQiaQjyQjj 

+ \m\AlAl + \m 2 CQ C£c£ + ±g 2 m (tft H ) (A% A® 

+ \glC*C*(Uy\ A \ B U*-e tl X t h t H l QiU*. (1) 

Here H is the Higgs doublet field, U(D) is the right-handed stop (sbottom) field, and Q is the 
third generation left-handed squark doublet field. The longitudinal components of the SU(2) and 
SU(3) gauge fields are denoted by A Q and C Q , respectively. The Latin (Greek) indices indicate 
SU(2) (SU(3)) components. As usual, the fields in eq. (|l|) are the static components of the scalar 
fields properly renormalized, the dimension of the fields in 3D is [GeV] 1 / 2 . Quartic couplings are 
of order g 2 (h 2 )T, having dimensions of [GeV]; here gi(h t ) denotes a gauge (top Yukawa) coupling. 
The zero-temperature trilinear coupling is X t = A sin j3 — fi cos (3 = A t sin f3. We work throughout in 
Landau gauge. In the following we have not included the correction to the quartic coupling between 
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the doublet Higgs field and the triplet scalar field A , or the corresponding correction for the SU(3) 
counterparts. We work throughout in the Landau gauge. We do not rewrite the expressions for the 
3D masses (rriA ,mc ) an d all of the quartic couplings which are not modified by the presence of 
the trilinear couplings, they are given in ref. 

2.1.1 Mass terms 

For the Higgs doublet we havef] 



9 2 L b Q -2 • 2 a L f 



m H:j = m H \l + -g w —-3h t smP m ^ 
X 3 , 1 



4 T ' 2 { 2 + 16^ + 1Q 9 ' 2 + 4 ^ Sin2 P + 4 (2/l * Sin2 P + 2X3 + K 



16tt 2 



6Am^ + 3(m 2 Q + m^h 1 sin 2 /3 + N c h 2 t X t 2 , (2) 



> 2 /2 \ 2 2 

where the Higgs mass parameter is denoted by run, and A = cos 2 2/3, A3 = ^f,A4 = —^f 

Similarly, for the third-generation squark mass terms we have 



< = <(l + 4, 2 ^) + T 2 (i, 2 + ^ + ^^in 2 /5 + ^ 2 ) 

~ I ^Q« + 2^^in 2 /5(m 2 ^ + m 2 3 ) + 2/ i 2 X 2 ), (3) 



- ^ {^qWq + 6A x m| + h\ml + h\ sin 2 /5m 2 , + h 2 t xf) , (4) 
where the soft SUSY-breaking mass for the third generation left-handed squark doublet is denoted 

2 2 

by mq, and X v = Ai = 

3 Throughout the paper we mostly neglect the hypercharge coupling g'. The only exception is in the contribution 
to the tree-level expression of the Higgs self-coupling A, as this latter quantity is fundamental in determining the 
strength of the phase transition. 
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2.1.2 Couplings 

We present here the expressions for the trilinear couplings in 3D, which were absent for the case 
of zero-squark mixing. The modifications of the scalar quartic couplings are suppressed at this 
stage for the values of the mixing parameters which we will consider, and we do not include the 
corrections from diagrams of those in fig. 1. However, at the next-stage of integration-out we will 
include them: 



X t h t = X t h t T^(l + 9 -g*-^--^h 2 t sm 2 p 



4 y 16vr 2 2 * H 16n 2 
+ ^^2) " X thtT 1/2 (\s + 2A 4 sin 2 (3 + 3h 2 sin 2 (3 + N c h 2 t 

~ 3*0 16^ (5) 
2.2 Second stage 

Another simplification of the effective theory can be obtained by integrating out the scalar fields, 
which are massive at the transition point. As we have seen, the static modes corresponding to the 
scalar fields Q, D, A a , C Q acquired thermal masses proportional to ~ g w ^T, as a consequence of the 
integration out of the non-zero Matsubara modes. The second stage proceeds in exactly the same 



way as in refs. |18| , |30| . We include the additional corrections arising from the couplings we have 
considered. 

2.2.1 Couplings 

The final expression for the tree-level 3D potential is given by 

V 3D = m 2 Hs H^H + \ H ,{H ] H) 2 + fn 2 Us U^U + X Us (U^U) 2 + %H^HU^U, (6) 
where the scalar couplings are now 

lb SirniAo o7TmQ 3 \ 
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+ (AD 2 ) + /if A? + /if AJ + i(^) 2 



+ 3X 2 /i t 2 (A^ + A* + h 



— 2 1.2 1 \ ( *> \ uL )—^- 1 

8 7rm Q3 (m Q3 + rmy 3 ) 2 



36 87rmc onm Ds \2 1 2 2 
L_ ( {h ?y - 2h? u g?" + 2h? U 9? 2 U + {9 U S1 D + 9 U J>) 2 + 2(£ D ) 2 



wm Q . 3 

x 2 Mh? u + 9?" + 9^)7^—1 ^ + X 2 t hhs: 



81 82 8irm Q3 (m Q , 3 + m H3 ) 2 8nm H3 {m Qz + m H3 ) 2 

- X t h$f(m Q3 ,m H3 ,m H3 ), 



7i = 73-^^(^ C7 + ^ 1 C/ + 3^ 2 l/ )(2A« + A2 + A: + ^) 

1 



and 



+ ^?(A? + £(AJ + /if))-4 



-2 



2 87rmQ 3 (m Q3 + m^ 



2 



+ x;h 2 t -(hr+gT + 3g 



2 V " 1 ' 1,81 82 87rm Q3 (m Q . 3 + m U:j ) 2 
3X 2 t h 2 t A Hs -^—- + AX 2 h 2 t A U3 - ' 



lTim Hz (m Q3 + m H3 ) 2 3 8nm U3 (m Q3 + m U3 ) 2 

X t h t ^—- ■ -- ■ -- ■ r- X t htf(mQ 3 ,m Ba ,m Ua ). 

47r (m Q , 3 + mu 3 ) (rriQ 3 + m Hz ) (m U3 + m H3 ) 



\ 1 2mi+m 2 + m 3 

/(m 1 ,m 2 ,m 3 ) 



^7r mi (mi + ms) 2 (mi + m2) 2 (m2 + ^3) 
2.2.2 Mass terms 



The one-loop contribution to the mass terms can be obtained directly as shown in ref. 0: 



W »> = < I 1 - 12 T ( m ;,' + m a) 3 j - ^'n,,, - ^(2A? + A^ + AS + h})m Q , 

~ 3h*X 2 t — L -, (11) 



- i^^^,-2^ M ^ +mm) . (12) 

In order to precisely fix the scales of the couplings that appear in the thermal polarizations of 
eqs. (2) and (3), one needs to perform a 2-loop evaluation of the effective potential. In addition, 
the mass parameters are renormalized in the 3D theory: 



m H 3 (ji) = m H 3 + hm H log (13) 



™u 3 (») = m U3 + ——f^ log (14) 
The expressions for the two-loop beta functions / 2mff , / 2m(l for the mass parameters have been 



given in ref. ||18|| . In appendix B we perform a two-loop calculation of the effective potential 
for the H (0-direction) and {/(^-direction) fields, including the effects of mixing with the third 
generation left-handed squark doublet^. We incorporate all of the corrections to the 3D couplings 
obtained in the previous sections, so as to determine the exact values of A# 3 and Ajj 3 - Using the 
results of appendix B it can be checked that the trilinear couplings do not produce scale dependent 
logarithmic contributions to the mass terms. In other words when both stop fields are light, the 
residual dependence on the mixing term in the divergent parts of the contributions to the mass terms 
appears only through the couplings which enter the expressions for the beta functions f2m H , f2 mu ■ 
These couplings are given by the one-loop expressions in equations (0) -® • We will analyse the 
effect of including these corrections on the critical temperatures £]. 

4 Some of the results given previously in the literature are corrected. 

5 The necessary modifications of fimn ? fimu m the case of non-zero mixing for large values of tuq appear in ref. 
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3 Results 



With the results of the previous sections and those in appendices B and C we can analyse the 
phase transition for values of mg = 300 GeV. In fig. |2| we show the critical temperatures for the 
transitions in the (solid)- and x (dotted)-directions as a function of the lightest stop mass m 4 ~ 2 , for 
tan/3 = 5 and for values of the parameter A t = = 100,200 GeV. We find that, for mg ~ 300 
GeV, there still is a region in which a two-stage phase transition can occur. This region is to the left 
of the crossing point of the curves. As shown in figs. § and |3], the range of values of the lighest stop 
mass for which the phase transition is sufficiently strong and a two-stage phase transition occurs 
decreases as X t increases. Thus the experimental bound on the stop mass puts an upper bound 
on the value of the mixing parameter. The figures display the effects produced by the non-zero 
trilinear couplings as described in the introduction. Keeping all other parameters fixed, lowering 
the stop mass will decrease the critical temperature in the 0-direction, slightly reduce the Higgs 
mass and strengthens the transition. The light-stop mechanism of enhancing the phase transition 
also eliminates the restriction of having low values of tan /3. 

The allowed region in parameter space is shown in fig. |], as a function of the Higgs and stop 
masses, for A t = 200 GeV, the region on the left of the solid line indicates when a sufficiently 
strong first-order phase transition occurs. The dotted line gives the condition for absolute stability 
of the physical vacuum. As explained in appendix C, to the left of this line the colour-breaking 
minimum is lower than the physical one at zero temperature. The dashed line is obtained when the 
critical temperatures of the transitions in the 0- and x-directions are the same. A two-stage phase 
transition occurs to the left of the dashed line. If the second transition from the CCB minimum 



to the electroweak minimum does not occur, as indicated by the analysis of ref. pi] , there are 
stronger constraints on the allowed region. For each value of the mixing parameter A t the available 
region is restricted to the band within the dashed and the solid lines. Note that, for low values of 
irtQ and non-zero squark mixing, there is no cross-over between the dashed and dotted lines. The 
end-points of the lines correspond to the maximum value of the Higgs mass which is reached by the 
effect of the zero-temperature radiative corrections for a given value of ijiq and A t . In fig. [5] we 
display the full allowed region in parameter space coming from the requirement of having a strong 



p0| |. The results of appendix B are no longer applicable as the high temperature expansion expressions for the D 
functions are not valid for the third generation left-handed squark doublet, see below. 
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enough phase transition. The solid vertical line is determined t varying tan/3 for A t = and the 
solid diagonal line corresponds to the variation of A t for tan/3 = 12. The region to the left and 
below these solid lines provides a sufficiently strong phase transition for values of 2 < tan (5 < 12 
and < A t < 280 GeV. The dashed line shown defines when the critical temperatures in the </>- and 
X-directions are equal for the same variations of tan/? and A t . For this value of uiq the maximum 
allowed masses are < 105GeV and m^ 2 < 170GeV. 

We now turn to the analysis of the phase transition for itlq = 1 TeV. The expressions for the 3D 
parameters are given in appendix A. The one-loop calculation is exact, the two-loop contributions 
are only estimated, using the corresponding two-loop beta functions for the mass parameters with 
non-zero mixing at large tuq given in ref. ||32|| . Here the scales A# 3 and Ajj 3 have been fixed to 



agree with the numerical results obtained in the case with zero- mixing in ref. |[30|| . We point out 
that the exact contributions from the two-loop sunset and figure eight diagrams when the high 
temperature expansion is no longer valid are not known []. These contributions are necessary for 
terms involving the third generation left-handed squark doublet for large values of tuq. Figure ^ 
shows that the phase-diagram structure is maintained for itlq = 1 TeV for several values of A t . 



Again, a possible two-stage phase transition persists, in agreement with previous results [f2(], |33 
We remark that the qualitative effects of increasing A t remain unchanged with respect to the case 
of mQ = 300GeV. In figure [7] the solid line determines the constant ratio |- = 1, as a function 
of the lightest stop mass for tan (3 = 5, varying the mixing parameter in the range < A t < 650 
GeV. The dashed line is obtained when the critical temperatures in the </>- and x~ directions are the 
same. The dotted line gives the zero-temperature condition for absolute stability of the physical 
minimum. Imposing the constraint which eliminatesthe region where tunneling into a CCB minima 
at finite temperature is possible reduces the allowed region to the band between the dashed and 
solid line. The experimental constraints on the stop mass m t ~ 2 > 80 GeV, will also restrict the value 
of the mixing parameter. Finally, in fig. [8] we present the contours of = 1 in the m^-m^ plane 
varying tan j3 for A t = 0, 100, 200, 300, 400, 600. The regions of possible two-stage phase transitions 
defined by each set of values of A t and tan (3 are not shown. The maximum allowed Higgs mass is 
now rrih < 120 GeV. 



3 We hope to return to this point in the future. 
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4 Conclusions 



We have performed a full two-loop dimensional reduction of 4D MSSM parameters to the 3D 
couplings and masses of the effective theory, including the effects from the trilinear terms in the 
potential. The preservation of the baryon asymmetry can be ensured tuning the stop mass in order 
to compensate the negative effects of the mixing term on the strength of the phase transition. We 
conclude that the direct stop searches restrict the amount of mixing in the stop sector. The allowed 
range of masses is m h < 107 GeV and < 170 GeV for rrig = 300 GeV. Our calculation allows us 
to determine the critical temperatures of the transitions precisely. We find that the phase-diagram 
still allows a possible two-stage phase transition for a small range of values of m^ 2 and any value of 
A t for values of tuq for which the high-temperature expansion can be applied. This range of values 
is shifted for different values of the stop mixing parameter. If the second phase transition does 
not occur there is a further restriction of the allowed regions of parameter space for electroweak 
baryogenesis for each value of the mixing parameter and tan f3. At large values of txiq the qualitative 
dependence on the parameters in the theory remains unchanged. 



A One-loop contributions from trilinear terms and large- 



The first diagram that contributes to the two-point functions corresponds to fig. la [34, |3"5| : 




A{m) + h{a) 



(15) 



where 




(16) 



and a 



Y- The expression that is valid for a high-temperature expansion is 




(17) 
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where L b = 2 log ^7 pa 2 log 7 Q ^ 5T , Lf = L b + 4 log 2. Here /I is the mass scale defined by the 
modified minimal subtraction (MS) scheme. For a low-temperature expansion that is m > T 
only the zero-temperature contribution remains as the temperature-dependent part is suppressed. 
Trilinear coupling induce an additional contribution to the two-point function shown in fig. lb: 

I 2 {m x , m 2 ) = B (p, m 1 ,m 2 ) + ^h( a , b), (18) 
where B Q (p, mi,m 2 ) is the usual Veltman-Passarino scalar function and in particular, 

B o (0,m 1: m 2 ) = -— A + 1 + In ^ + 1 In ;! 



167r 2 V m 2 \m\ — m\) m 2 , 

2 1 Ji 



/ 2 (a,6) = — ^ dx( - j r j (19) 

a 2 -b 2 Jo \(x 2 + a 2 )5 e (*"+a*)3 _ i (x 2 + b 2 )2 e {x^)h _ X J 



where a = mi/T, b = m 2 /T. 

For a high-temperature expansion this becomes 



h(mi, m 2 ) = JL(1 + L 6 ) + O(^) + O(^). (20) 
For mi < T and m 2 > 2ttT is 

I 2 (m 1 ,m 2 ) = B (p,m 1 ,m 2 ) + -\fi{m 2 ) (21) 
and for the case in which mi = m 2 ^> T this reduces to the zero-temperature contribution 

1 u 2 

I 2 (m,m)=B o {0,m,m) = -—ln^- (22) 

1D7T rrr 

plus wave-function renormalization effects. We remind the reader that this can be obtained from 
7 2 (m, m) = —q^t I(m). Similarly, we can obtain the corresponding results for other diagrams, in 
fig. 1 taking the appropiate derivatives. We now give only the low-temperature expressions valid 
for mi T, the corresponding high-temperature expressions can be found in the literature: 

I | Tfl 2 

7 3 (mi, m 2 , m 2 ) = C o (0, m 1 ,m 2 , m 2 ) + zr^— In (23) 

107T 771-^ jj,rp 
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I 3 (mi, mi, m 2 ) = C o (0, mi, m 1; m 2 ) + ^r/i(m 2 ) 

mf 



I A (m l} m 1} m 2 , m 2 ) = -D o (0, mi, mi, m 2 , m 2 ) + 



1 1 ^ m\ 
mf 167r 2 /iy 



where the functions C and D a are the scalar one-loop functions of Veltman-Passarino. 
The 3D couplings in equation ^ for the case of large rriQ are given by: 



6fc? sin 2 /? 



2^16tt 2 



16tt 2 / 



16tt 2 



L b (j^ 4 +12A 2 + i/i 4 sin 4 /3) 3[ \i, -\ A :! A , -f A^ cos J ;+A 2 sin J / 



+ h 2 sin 2 (3(X 3 + A 4 sin 2 (3) + \h A t sin 4 (3) In + + 3/i 4 sin 4 

2 / m,^) 8 



+ 3X 2 /i 2 ( A 3 + A 4 sin 2 (3 + hf sin 2 (3) 



m 



Q 



777 

+ 3X?htsm 2 {3(--^(l-ln^-))-^X?hti 



3 
2" 



1 



4 (2 + ln4))-^ 2 4 



Al/3 " l T ( 1 + 8 ^I&)"rD^ 



+ f + /> 4 sin 4 /?) In 4 + ^ + 2X 2 /, 2 (/, 2 - k 2 ) -L 
o / mg 18 3 



2X 2 /i 4 sin 2 /?( o-(l - In 



m} 



■(2 -In-*)) 



h, 

3 'm^ 



7 3 



/i 2 sin 2 /3T 1 + T # 



167T 2 



4 ^l 67r 2 ^ sm l 67 r 2+4 ^167r 2 J 
L 6 f-/t 2 sin 2 /5c/ 2 + 2/i 4 sin 4 /5 + 6Xh 2 t sin 2 ^ 



( h 2 (2X 3 + A 4 + /i 2 sin 2 £ ) In \ + X 2 t h 2 t (2(A 3 + -A 4 sin 2 /? + h 2 sin 2 /3 + h 2 t )\ + (6A 

mn \ 2 TTln 



+ 2/> 2 sin 2 fS){-\{! - In ^)) - X 2 h 2 t (-^(2 + In ^)) - ^ sin 2 /5(^ + 4 
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and the 3D masses are given by 

Q 

_ J±- f 6 Am 2 + 3m 2 )A 2 sin 2 fll - f_^L + < (fa A 4. 1^ 

+ 1„ £l + 1) - - J-*,™* + ^/ 2m „ log ^, (29) 



< - m K 1+ ^ife) +r2 (^ + ^ + ^ sm ^)- 2 '#( to £5 +1 ) 

u 2 m 2 rj f, at \\ ,,X 2 T 2 1 1 A™, , , 

+ '»£j +i +4( ln £| + v)- 2ft '4^-^ (30) 

where \it = 47re~T. 



B Two-loop contributions with non-zero squark mixing 

The strategy we employ follows that of refs. J/], The idea is to use the 4D two-loop effective 
potential in order to fix the scales in the 3D theory, and to use the 3D effective potential expressions 
for the Higgs and stop fields given in ref. |18| to analyse the phase transition. We calculate 
the unresummed two-loop effective potential in order to include all 4D corrections to the mass 
parameters; resummation is automatically included in the calculation of the two-loop effective 
potential in the 3D theory. We must also include the contributions to the two-loop effective potential 
of the static modes, which have been integrated out at the second stage (includes the effects of 
resummation of the heavy fields). 

There are several effects that must be considered in order to obtain all of the contributions 
(constant and logarithmic) to the mass parameters. From the 4D effective potential, one finds the 
two-loop contributions from the gauge bosons, Higgs, right-handed stop, left-handed squark doublet, 
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right-handed sbottom, and top quark. The main difference is that the D functions appearing below 
correspond to the unresummed expressions. Additionally we must include the effects arising at the 
second stage of reduction from the left-handed squark doublet, the right-handed sbottom, the scalar 
triplet and the scalar octetQ 

We now derive the effective potential at finite temperature using the background fields <ft an d 
X = tR a u a , where we have chosen the unit vector in colour space u a = (1, 0, 0). We first write the 
expressions in the shifted theory of the mass spectrum after the first stage of integration. 

B.l ^-direction 

The gauge-boson masses are 

m w,z = -^9l4> 2 - (31) 
The stop mass matrix elements are given by 

ml((f>) = m 2 Q + ml{4>) + \{\ - jj sin 2 9 w ){g 2 w + g' 2 )<p 2 cos 2/? (32) 
™U<I>) = m l + ™ 2 t (<f>) + \i\ sm 2 w ){g 2 w + g' 2 )<p 2 cos 2/5 (33) 



m i LR ^ = sin/? -ficoB/3)<f> = ( 34 ) 

the corresponding eigenvalues are and m t - 2 , and the eigenstates are given by 



ti = cos a t tL + sin a t tji 

t2 = — sin a t ti + cos a t tR, (35) 

2m? 

where sin2«t = , 2 , . Below we will use the abbreviations c t = cosa 4 , s t = sinc^, etc. 
_ h m h' 

Neglecting mixing effects in the sbottom sector, we have 



m\ L (<f>) = m\ + \{-\ + ~ sin 2 6 W ) (g 2 w + g' 2 )^ 2 cos 2f3 I :$(> ) 



The expression given for the 4D effective potential would correspond to the usual resummed two-loop 4D effective 
potential if we used the resummed expressions in the D functions appearing below. 
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ml R {4>) =m^ + i(-isin 2 M(^ + ^ 2 )0 2 cos2/5. (37) 
For the Higgs sector, the Goldstone bosons and Higgs masses are 

2 2 , \ i 2 

m * = m H + X(f> , 

m 2 = m | + 3A0 2 . (38) 

The additional corrections that arise from supersymmetric particles can be calculated using the 
two- loop unresummed potential. Our notation for the D-functions corresponds to that of ref. 0. 
The contributions from the two-loop graphs containing supersymmetric particles are given below. 
For the (^-direction, we can drop the colour index of the squark masses: 

g 2 

(SSV) = — ^-N c [ctD ssv (mi , m iv m w ) + s^D ssv (m i2 ,m i2 ,m w ) + 2s t c 2 D ssv (m il ,m i2 ,m w ) 

o 

+ D SS v(m- bL ,m- bL ,m w ) + A(c 2 Dssv{rn- t ^ m b L ^ m w) + s 2 t D ssv (m h ,m bL ,m w ))] 
- ^(N 2 - lj^K-^-.O) + D ssv (m h ,m h , 0) 

+ D ssv (m {2 ,m i2 ,0) + D ssv (m~ bR ,mi R ,a)], (39) 



+ 



N 

(SSS) = -f 



h 2 sin 2 (3 + ^c 2 cos 2/3)0+ V2X t h t c t s t 



+ 
+ 2 



/i 2 sin 2 /3 + -j-s 2 t cos2[3j(j) - V2X t h t c t s t D S5S (m t - 2 , m t - 2 , m h ) 
(9 2 \ 1 I 2 



+ (^f cos 2/30] D sss (m~ bL ,m~ bL ,m h ) 



+ h 2 t X 2 D SS s (mfc , m fa , m T ) + 



/i 2 sin 2 /3 + ^ cos 2/3^ q + V2h t X t s t 



1 2 



/i 2 sin 2 /3 + ^ cos 2/3 ) s t - V2h t X t c t 



Dsss(m i2 ' m b L ' m K 



Dsss{m il ,m b m n ) 



(40) 
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(SV) = -\g 2 s {N 2 c - l)[D sv {m h ^) + %K,0) + D sv (m i2 , 0) + L> sy (m- 6fl , 0)] 
3 

- ^9 2 w N c [c 2 t D sv (m h) m w ) + s 2 D sv (m h , m w) + D sv (m bL , m w )\, (41) 



(55) = ~^N C {N C + l)[4Z>ss(TO fl ,m fl ) + D ss (m~ bR ,m~ bR ) + c 2 2t D S s(m i2 ,m i2 ) + D S s(m bL , mgj] 

+ jN 2 (ct + st + lQc 2 s 2 t )D ss {m h ,m h ) - |iV c (^ + sf - 2 -s 2 c 2 t )D ss {m iv m h ) 

+ ^iv c (2 - AQ[c 2 L>s S (m fl ,TO S J + s 2 £>ssK~ 2 , m 6 J] 

+ h 2 N c [(c 4 t + s*)D ss {mt v m~ t2 ) + s 2 t D S s{m bL ,m~ ti ) 

+ c 2 D ss {m- bL ,mt 2 ) + {N c + l)c 2 t s 2 t [D S s{m- tv m h ) + D ss {m~ t2 , m h )\ - 2N c s 2 c 2 Dss{m i2 ,m h )} 



+ 



)iV c (iV c + l)[c A t D ss {m k ,m k ) + 2c 2 S 2 D S5 K 1 , m t ~J 



+ s^Dss {m- t2 ,m {2 ) +D S s(rn bL ,m~ bL ] 



+ N, 
+ 



1 1 

-h 2 sin 2 (3 + -gl cos 2/3c 2 ) [D S s{m k > Wto) + D ssi.m- tl , 



1 1 

-/i 2 sin 2 /3 + -# 2 cos2/3s 2 ) [D ss (m t - 2 , m ft ) + D ss (m~ t2 , m w )] 



1 



- ^c9 w cos2/3[Dss(m- bL ,m h ) - D ss (m~ bL ,m^) + 2c t D ss (m~ tl , m n ) 
+ N c h 2 sin 2 PlstDssimi^mn) + c 2 D ss {m- t2 ,m 7T ) + D S5 (m bi , m^)]. 



2s 2 L> 5S (m 4 - 2 ,m w )] 



(42) 



There are no additional finite contributions apart from those given in ref. [g(J, from counterterms 
when non-zero mixing is included in the stop sector. One can check that when all contributions are 
added there is no dependence on the mixing angle in the divergent part of the potential. 



B.2 ^"direction 

The gauge-boson masses are 



2„,2 



m G = ~ A 9sX 



III 



G 



4 2 

- 3 m G . 



(43) 
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For the Higgs sector, the Goldstone bosons and Higgs masses are 



2 2 2 i \ 2 i u2 ■ 2 n't* 

ml = ml + 3\ X 2 + h 2 t sm*(3^. (44) 



The masses of the rest of the scalars contributing to the effective potential are given by 



m 2 hl =m 2 H + h 2 t sm 2 p^, (45) 
<i = ~^ X *X ( 47 ) 



m 2 h2 =m 2 H + h 2 t sin 2 (48) 



,2 _ _2 / r,2 9s\X 2 



<=-Q+(5-fJy> (49) 



so the mixing angles are related by 



with eigenstates given by 



sin qi = — sin q 2 = s g (51) 



#1 


= Cg/ll - Sgb Ll 


92 


= S q hl + C q b L1 


#2 


= C q k 2 + Sg^Ll 


92 


= S q h 2 + C q ti 



(52) 
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TO? 

*L 9 



(53) 



TO 



1>L 



2.3 



(54) 



2 2 , /^AX 2 
TOr = TO n + , 

bR i D \ 6 / 2 



(55) 



2 2 ^ 2 \X 2 

TOr = TO n — . 

6 «2,3 D V 3 / 2 



(56) 



In the expressions below, we also include the contributions from the Higgs doublet. The two-loop 
unresummed effective potential in the x-direction is given by the following contributions 8 : 



g 2 

(SSV) = — ^-[D S sv{m qi ,m qi ,Qi) + D S sv(m Hl ,m Hl ,0) + D S sv{m L2 ^ m t L2 

o 

- D ssv (to,-, „ , mr, , 0) + D ssv (m q2 , m q2 , 0) + D ssv (m H2 , m H2 , 0) 



,0) 



+ Dssv(m iL3 ,in~ tL3 , 

mi. 3 ,0) + D SS v(m~ bL3 , m- bL3 , 0) 



+ A(D S sv{mq l ,m h ^) + D S sv(m Hl ,m H2 ,0) 



+ L> 55 y (TO t - i2 , TOj l2 , 0) + L> 5S y (m iL3 , TOj L3 , 



0))] 

(iV c - l)Z?s S y(m w ,m w ,m G ) + (iV c - l)D 5S y(m w , m u , m G ) 

N c - 1 ^ ._ 1 _ . 

+ — — Ussvim^, m u , m G ) + —D S sv\m w , m u , m G ) + 

+ 2(N C - l)(c 2 q D S sv{m qi ,m iL2 ,m G ) + s 2 q D S sv(m H2 , m iL2 , m G )) 

^^(4 D ssv(m qi ,m qi ,m G ) + s q D S sv(m H2 ,m H2 ,m G ) 



1 

m G ) + N c(N c - 2)^55^(^0;, m w 



,0) 



+ 



^4 c l D ssv(m qi ,m H2 , m G )) + —D S sv(m iL2 ,m iL2 ,m G ) 

< m t L2 ^) 



+ N c (Nc-2)£>ssvK M , 



+ 2(iV c - l)(c 2 q D S sv(m q2 ,m bL2 ,m G ) + s 2 q D ssv (m Hi: m bL2 ,m G )) 

= ni( tj , mj i2 = TOg , mi R2 = m l R3 m tnc x-direction, wc just multiply by a factor of 2 the contributions 



8 As nif — mf ,mr =mr ,mr = mr m th< 
from these fields in some of the following expressions. 
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+ ^ 



r c — 1 / \ 
^ — \c A q Dssv{mq 2 ,mq 2 ,m G ) + s q D S sv(m Hl ,m Hl ,m G ) + 2s 2 q c 2 q D S sv(rnq 2 ,m Hl ,m G )J 

+ -j^D SS v{m~ bL2 ,m~ bL2 ,m G ) + N C (N C - 2)D SS v(m~ bL2 , m bL2 , 0) 

iV c — 1 

2(iV c - lj^K^.^^c) + ^- D S sv(m bi 



+ 



+ ^ D 55v(m- bfl2 , m- bR2 ,m G ) + N C (N C - 2)D SS v(m bR2 ,m bR2 , 0) . 



(57) 



(SV) = -^[2N c (N c -2)D sv (m iV ,0) + (N c -l)[3D sv (m iV ,m G ) + D sv (m u ,m G )] 



+ —[(Nc+^Dsvim^rnG) + {N c - l)D sv (m u ,m G )] 

+ 2N C {N C - 2)D sv (m iL2 ,0) + (N c - l)[2D sv (m iL2 ,m G ) + 2(c 2 q D sv (m qi ,m G ) 

+ s 2 q Dsv(m H2 ,m G ))} + ■^-[2D sv {'m iL2 ,m G ) + 2(c 2 q D sv (m qi ,m G ) + s 2 q D sv (m H2 ,m G ))} 

+ 2N C (N C - 2)D sv (m bL2 ,0) + (N c - l)[2D sv (m bL2 ,m G ) + 2(c 2 q D sv (m h ,m G ) 

+ s 2 q D S v(m Hl ,m G ))} + — [2D sv (m bL2 ,m G ) + 2(c 2 q D S v(mq 2 , m G ) + s 2 q D sv (m Hl ,m G ))} 

+ 2N C {N C - 2)D sv (m bR2 ,0) + (N c - l)[2D sv {m bR2 ,m G ) + 2D sv (m bRi ,m G )} 
1 

- ^[%K,0) + D sv (m iL2 ,0) + D sv (m iL3 ,0) 

+ D sv (m q2 ,0) + D sv (m~ bL2 ,0) + %(mj M ,0) + D sv (m Hl ,0) + D sv (m H2 ,0)}, 



— [2D SV (m bR2 , m G ) + 2D SV (m~ bm , m G )] 



(58) 



(SSS) = -\jjX [3D S ss(m u , m u , m u ) + (2N C - l)D S ss(m u , m w , m w )] 

If 1 \ 2 

- - [h 2 sin 2 /3cjx + (h 2 - ^g 2 s )s 2 q x + V2X t h t c q s q j [D sss (m u , m Hl ,m Hl ) 

+ D sss (m u ,m H2 ,m H2 )] 

If 1 \ 2 

- 2 ( h t sin2 P s v x + ( h * ~ ^9s)c 2 q X ~ V2X t h t c q s q ) [D SS s(m u , m qi ,m qi ) 

+ D sss (m u ,m q2 ,m q2 )] 
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/ 1 y \2 

- 2 (^t - 2 9s " ) ^ Cq ~~ s i Xtht ) [-°S55(^a;, m qi ,m iL2 ) + D SS s(m w , m q2 , m bh2 ] 

- 2 [0t - + c g X t ^] 2 [D 5 55(m w ,m^ 2 ,m t - i2 )) + D S ss("i«,"iHi,"i 6j . 2 )] 

.AT 1 

- (-/i 2 sin 2 /3s g CgX + c 2g ^= + s g c g (/i 2 - -# 2 )x) 2 £>sss(m«, m H2 , m qi ) 

X 1 

- (h 2 sin 2 (3s q c q x + c 2g -/= - s ? c ? (/i 2 - -g 2 )x) 2 D sss (m u ,m Hl ,m q2 ) 

- [c2g y= h t ] 2 (D S ss {mu , m Hl , m q2 ) + D SS s (m^ , m H2 ,m qi )) 

- ^(2(^) 2 Dsss(m u ,m iL2 , mi J + 2(^) 2 D sss (m u ,m bL2 ,m bL2 ) 



,9 



+ 2d) 2 D sss (m u ,m- bR2 ,m- b J), (59) 
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= 2A c q (D S s(rn hl ,m hl ) + D S s(rn h2 ,m h2 ) + D S s(rn hl ,m h2 )) 

+ s q (D S s(m qi ,m qi ) + D ss (m q2 ,m q2 ) + D S s{m qi ,m q2 )) 

+ c 2 q s 2 q (2D S sijn hl ,m q2 ) + 2D ss {m h2 ,m qi ) + D ss (m hl ,m qi ) + D ss (m h2 ,m q2 )) 
( g 2 g 2 \ 

+ i-f + )[2c q D S s(m qi ,m qi ) + 2s q D S s(m h2 , m h2 ) + As 2 q c 2 q D S s{m h2 ,m qi ) 

+ 4(sjD ss (m fc2 , m fL2 ) + c 2 q D ss (m qi ,m iL2 )) + 6^55(7^^, m t ~ M ) 

+ 2c q D S s(m q2 ,m q2 ) + 2s q D S s{m hl ,m hl ) + 4s 2 q c 2 q D S s(m hl ,m q2 ) 

+ A(s 2 q Dss{rn hl1 m iL2 ) + c 2 q D ss (m q2 , m iL2 )) + 6D ss (m bL2 ,m bL2 )} 
1 3 2 

+ -(-T-)[SD S s(m u ,m u ) + l(XDss(m u , mj + 35Ass(?tv,, m^) 
4 o 

+ 8£) S sK m ,m gm ) + 2AD ss (m bR2 ,m bR2 ) + 16^^,™^)] 

+ _ ^dl^qDssimq^mu) + c 2 q D ss (m qi , mj + s 2 q D S s(m h2 ,m u ) + s 2 q D S s{m h2 ,mS) 

+ 2D ss { , mi L2 i m u) + 2D ss (m iL3 , rn^)] 

+ ^^[^Ass^^m^) + 5c 2 D 5 s{mu,m qi ) + s 2 q D S s(m u ,m h2 ) 
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+ 5sj£>ss(m w , m h2 ) + D ss (m u , m~ tL2 ) 

+ 5D SS (m w , 77if M ) + D ss (m u , m iL J + 5D ss (m u , m^J] 

+ - ^s)[^-Dss(^, 2 , m u ) + (? q D ss {m^ m u ) 

+ s 2 q D ss (m hl ,m u ) + s 2 q Dssijn hl ,m w ) + 2D ss {mi L2 ,m U} ) + 2D ss (m iL3 , mj] 

+ ^# 2 [ c gAss(™ u , wiga) + 5cjD 55 (m w , m g ~ 2 ) + s 2 q D ss (m u , m hl ) 

+ 5s^55(m w , m hl ) + D S s(m u , m~ bL2 ) 

+ bDssimuj, m~ bL2 ) + D S s(m u , m bh3 ) + 5-Dss(m w , ™g M )] 

+ (/i 2 sin 2 /3 + ^ 2 cos 2/3) 2c 2 s 2 [ J D 5 s(m fel ,m ftl ) 

+ D ss (m h2 , m h2 ) + D ss (m hl , ) + D ss ijn qi , m q2 ) + D ss ijn qi , m qi ) 

+ D ss (m q2 ,m q2 )} + c 2 q [D ss (m iL2 ,m h2 ) + D ss (m iL3 ,m h2 ) 

+ D S s(m~ bL2 ,m hl ) + D ss (m bLs , m hl )} + s 2 q [D ss (m iL2 , m qi ) + D ss (m iL3 , m qi ) 

+ D S s(m bL2 ,mq 2 ) + D ss (m bL3 , rriq 2 )] + ^Jl^fm^, m fe2 ) + D ss (m q2 , m hl )\ 

+ (^w( c J + 4 - ( c t + s i)\dl cos 2 /? - 2s 2 c 2 (-^ 2 sin 2 /3 + h 2 sin 2 [£>ss(m 5l , m ftl ) 



g 2 

+ /^(m^m^)] - ysJc^Ds^m^m^) + D S s{m qi ,m q2 )\ 
~ {j9icos2(3^[c 2 q {D S s{m- bL2 ,m h2 ) + D ss (m bL3 ,m h2 ) 

+ D ss {m- tL2) m hl ) + D ss (m- tL3 , m hl )) + s 2 q (D ss (m bL2 , m^) + D ss (m bL3 , m^) 
+ D ss {m~ tL2 , m g ~ 2 ) + D ss {m~ tL3 , m q2 ) )] 

+ sin2 Pi c t D ss( m u, m hl ) + c 2 q D ss {m u , m h2 ) + s 2 q D ss (m u , m qi ) + s 2 q D ss (m u , m q2 ) 



+ 



5(c 2 D 5 5(m w ,m /il ) + c 2 q D S s{mu,m h2 ) + s 2 q D S s(rriu, m^) + s 2 L> 5 s(m w , TOg 2 ))] 



+ ^(2-7V c )[c^ 55 (m ? - 1 ,m ? - 2 ) 



+ s q D S s{m hl ,m h2 ) + s 2 q c 2 q Dss(m h2 ,m q2 ) + s 2 q c 2 q D S s(m hl ,m qi ) 
+ D S s(m iL2 ,m bL2 ) + A?s(mf iS , m^)] 

+ 2 9l\ c \ D ss (rriq^m^) + s A q D S s{m hl ,m h2 ) + s 2 q c 2 q (D ss {m qi ,m hl ) 
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+ D ss (m q2 , m h2 )) + D ss (m iL2 , m~ h J + D ss (m~ bL3 , m h3 )\ 

- ^g 2 s [Cq D ss{rn qi ,m q2 ) + s q D ss (m hl ,m h2 ) + s 2 q c 2 q (D ss {m qi , m hl ) 
+ D ss (mq 2 ,m h2 )) + clD S s{mq^mi L2 ) + s 2 q D S s{™.h 2 , m~ bL2 ) 

+ c 2 q D ss {m qi ,m~ bL3 ) + sjD 5S (m h2 ,m^ 3 ) 

+ c 2 q D ss (m iL2 ,m q2 ) + s 2 q D ss (m iL2 ,m hl ) + D ss (m iL2 ,m bL2 ) + D ss (m~ tL2 ,m b J 

+ c2 q D ss(m iLS ,m q2 ) + s 2 q D ss (m iL3 , m hl ) + D ss (m iL3 ,m bL2 ) 

+ D 5 5K i3 ,m 5M ) 

- ^9 2 s [c 2 q D ss (m qi , m bm ) + s 2 q D ss (m h2 ,m bRi ) + D ss (m iL2 ,m bR2 ) 
+ D S s(m iL3 ,m bm ) + c 2 g D ss (mg 2 ,m bm ) + s 2 q D S s(m hl ,m bm ) 

+ D ss (m bL2 ,m bR2 ) + D ss (m bL3 ,m bR3 )} 

+ -Q9 2 s{c 2 q D ss {m qi ,m bm ) + s^D^m^, m bm ) 

+ <? q D S s(™.h, m bR2 ) + s^55(^ 2 , m Bj J + (? q D ss {m^, m bR3 ) + sjD 55 (m ft2 , m Bj J 

+ ^55 {rn~ tL2 , m Bm ) + D ss {m~ tL2 , m bR2 ) + D ss {m~ tL2 , m bR3 ) 

+ D ss (m iL3 , m bRi ) + D ss K i3 , m bR2 ) 

+ ^55(^3,^3) 

+ c 2 q D S s{mq 2 , m bm ) + s^D^m^, m bm ) 

+ cjD 5S (mg 2 ,m^ 2 ) + s^D^m^m^) + c 2 q D ss (m q2 , m bR3 ) + sjD 5S (m ftl , m bR3 ) 

+ D ss (m~ bL2 ,m~ bm ) + D ss (m~ bL2 ,m~ bR2 ) + D ss (m bL2 ,m bR3 ) + D S s(m bL3 ,m bm ) 

+ D ss (m bL3 ,m bR2 ) + D ss (m bL3 ,m bR3 )} 

+ -g 2 [D S s{m u ,m bm ) + D 55 (m w ,m Bm ) + 2D ss {m ul , m bm ) + 2 j D 55 (to w , m Sj J] 

+ 5^ 55 (m w ,m Sm ) + 5D 55 (m w ,m^ 2 ) + hD S s{m w , m bR3 )\. (60) 
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B.3 Integration over the heavy scale 

The second part of the calculation arises, as noticed in the paper by Kajantie et al. 0: when the 
"heavy" particles have been integrated out their contributions to the 3D mass parameters should 
also be included, as they can substantially vary the value of the parameters A^ 3 ,A[/ 3 . In order 
to do this we must calculate the two-loop contributions to the effective potential in the 0- and 
X-directions from the heavy fields: Q, D, C a , A Q . In the following rotation to eigenstates the angle 
is temperature-dependent; however, if we verify that the eigenstates are always well separated close 
to the transition point, there is no ambiguity about which is the field that is being integrated out 
at the second stage. 

The masses in the shifted theory are now given by 




(61) 




(62) 




(63) 




(64) 




(65) 




(66) 



and the relevant mixing terms are: 




(A t sin (3 — fi cos 0)<j) = —=X t (f) 




(67) 



for the 0-direction, and 



mU<f>) = ±^X tX 
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for the x-direction. 

The expressions for the rest of the fields are given in [IB]. The two- loop contributions from 
the heavy scale are given below. We stress that the Z)-integrals in eqs. and ( [70| ) are just 3D 
integrals, our notation follows that of refs. |TE, [7| 9 . We do not write the contributions arising 
from the longitudinal components of the gauge fields A and C Q , since the only modification that 
is necessary is to substitute the gauge couplings g w ( s ) 9w(s) 3 , which are given in ref. |30 



B.3.1 ^-direction 



g 2 

yheavy^ = — ■j*-N c [c$D sgv (m il ,m il ,m w ) + StD SS v{mt 2 , m t 2 ' m w) + 2s 2 t c 2 D ssv (m~ ti ,m i2 ,m w ) 
?ssv(m~ bL ,m~ bL ,m w ) + A(c 2 D S sv{m { n m b L i m w) + s t D ssv(mt 2 i m b L i m w))] 



+ D t 



9 -f(N 2 c - l)[ctD ssv (m h ,m h ,0) + stD ssv (m h ,m i2 ,0) + 2cls 2 t D ssv (m h ,m i2 ,0) 
+ D S sv{m- bL ,m- bL ,0) + D SS v{m- bR ,m- bR ,0)] 

\h L t + A 3 + K)4)<t> + V2X t h t c t s t 



N. 



D 



+ 

+ 2 



\ 1 

(hf + A 3 + A*)s t 2 U - V2X t h t c t s t D sss (m- t2 , m i2 ,m h ) 



{h L t + A 3 + A s 4 )c t s t ^j4>+ -f=X t h t c 2 t 

2 



D 



^sssimt^m^rrih) 



[hz + A-l) <pD 
+ h 2 X 2 t D 



+ 



'sss(m bLl m bL ,m h ) 
h 2 X 2 t D S ss{m h , m h ,m % ) + - A^ + A^ c t + V2h t X t s t 
\f -At + Aljs t - V2h t X t c t l2 



D 



'sss(m- tl , m b^ m ^) 



D 



- \g 2 S3 {N 2 c - l)[4D sv (m k ,0) + s 2 D sv (m h , 0) + D sv {m h , 0) + D sv {m~ bR , 0)] 
3 

_ Qgl^clctDsvimt^mw) + +s 2 D sv (m i2 ,m w ) + D sv (m bL ,m w )] 
+ N C (N C +1) [\ Q: 4 + g^s 2 t c 2 t + gf 2 u c 2 s 2 ] D ss {m k , m {l ) 

+ N c (N c +l)[\ Qr s A t +gi u s 2 t c 2 + gf 2 u c 2 s 2 }D ss (m i2 ,m i2 ) + Xq 3 N c (N c + l)D ss (m~ bL ,m iL ) 
9 Our convention for the functions Dyvv, Dyvs is that of |is}| . 
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+ N?(2c^ t \ Q . A - 2c 2 s 2 gf l u + (of + st)g? 2 u )D ss (m k ,m i2 ) 

+ N c (2c 2 s 2 \ Q3 - 2c 2 s 2 gf 2 u + {c\ + s\)gf)D ss {m ix , m h ) 

+ \ D3 N C (N C + l)D ss (m- bR ,m- bR ) + (2A 1 )N C (2 - N c )[c 2 D ss (m- ti} mJ 

+ s 2 t D ss (m i2 ,m bL )} 

+ (N c g^ + N 2 g^)[c 2 D ss (m ii: mJ 

+ St D ss(m i2 ,m bL )] 

+ N e {gf + N c g^)[s 2 t D ss {m^mJ + c^K, mj] 

+ h? u N c [(c 4 t + st)D S s{m h ,m h ) + s 2 D ss (m bL , m~ tl ) 

+ c 2 D ss {m lL , m h ) + (N c + l)c 2 s 2 [D ss (m h , m~ tl ) + D ss (m~ t2 , m h )] - 2D ss (m i2 , m~ tl )] 

+ A 1 N c (N c + l)[c 4 t D ss (m il ,m il ) + 2c 2 s 2 D ss (m h ,m i2 ) 

+ s*D S s(m i2 ,m h ) + D S s(m~ bL , m~ bi )} 

+ 



(h L t + A 3 + A s A )s 2 j [D ss (m- t2 , m h ) + D ss (m h ,m n )} 
^N C (A 3 + A c 4 )[D S s(m~ bL ,m h ) - Dss(m bL ,m w ) + 2c 2 D S s{.m- t ^ m ^) + 2s 2 Ass(m t ~ 2 , m n)} 



+ N c hf[Dss(rn bL ,m v )] 

+ (<7?^ + 3gf 2 D )N c {c 2 D S s(m- tl , mj + s 2 D S s(m i2 ,mJ + D ss (m bL ,mJ 
+ (<C + 3^ 2 D )7V c [ St 2 L> 55 (m fl , mj + c 2 D ss (m i2 , m-J}. 
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B.3.2 ^-direction 



yheavy^ = — m.^A Dssv ^ m ^ ? m ^ ? q) + Sq D S sv{m Hl , m Hl , 0) + 2c 2 q s 2 q D SS v{m qi , m H2 , 0) 

+ £W(m t ~ i2 , m- th2 , 0) + D S sv{m~ tL3 , m im , 0) 

+ CgDssv(mq 2 , mg 2 , 0) + s q D S sv(m H2 ,m H2 , 0) + 2c 2 s 2 ,LW(ra ( j 2 , 0) 

+ r> 55y (m gj . a , m^ 2 , 0) + r> 55y (m^ 3 , m^ 3 , 0) 

+ ^(CgDssvimq^m^^) + s q Dssv(m Hl ,m H2 ,0) + 2c 2 q s 2 q D SS v(m qi , 

+ 2cVr 



{c q ussv {jn^ , m h , u; -+- s q u S sv \jn Hl , 

c 2 q s 2 q D S sv(mq 2 ,m H2 , 0) + D S sv(mt L2 ,m bL2 , 0) + D S sv(m iL3 , 



.q^mn^O) 
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- 9 2 S3 \[ 2 ( N c - l)((^D ssv (m 4l ,mi L2 ,m G ) + s 2 q D S sv(m H2 ,m iL2 ,m G )) 
c q D S sv(m gi ,m qi ,m G ) + s q D ssv (m H2 ,m H2 ,m G ) 



N c -1 

+ 



2s 2 q c 2 q D S sv{m qi ,m H2 ,m G )) + —D S sv(m~ tL2 , m t L2 ,™G) 



+ N c (N c -2)D SS v(m iL2 ,m iL2 ,0) 

+ 2(N C - l)(c 2 q D S sv(m q2 ,m- bL2 ,m G ) + s 2 q D S sv(m Hl ,m~ bL2 ,m G )) 

+ N ° N 1 (c q D S sv(mq 2 , m q2 , m G ) + s q D S sv(rn Hl ,m Hl ,m G ) + 2s 2 q c 2 q D S sv(m q2 ,m Hl ,rn G ) 

+ -^D SS v{rn~ bL2 ,m~ bL2 ,m G ) + N c (N c -2)D S sv(™b M > m b M >°) 

+ 2{N C - l)D S sv(m~ bRi ,m~ bR2 ,m G ) + N ° 1 D S sv(m bm ,m~ bm ,m G ) 

+ jj- D ssv(m~ bR2 , m bR2 ,m G ) + N C (N C - 2)D SS v(m bR2 , m bR2 , 0)] 

+ + 9® 2 U )s 2 q X + V2X t h t c q s q ) 2 [D S ss(m u , m Hl ,m Hl ) + D SS s(m u , m H2 , m H2 )\ 

~ \(( h ? U + 9^ + 9? 2 U )4x - V2X t h t c q s q ) 2 [D S ss(m u , m qi ,m qi ) 

+ D S ss(m u ,m q2 ,m q2 )} 

- 2[(h? U + g®V)— Cq - s q X t ht] 2 [Dsss(rn u , m qi ,m iL2 ) + Dsssim^, m q2 , m b J 

- 2[(hf u + gf^^Sg + CgXthtflDsssim^.mH^mtJ) + Dgssim^, m Hl , m s J] 

- (c 29 -^= + s q Cg(h? U + gf? + gf 2 U ) X ) 2 D sss (m u , m H2 ,m qi ) 
~ (c 2q ^ - s q c q (hf u + g?u + gf 2 U )x) 2 D S ss(m u , m Hl ,m q2 ) 

- [c2 q ~^ h t\ 2 ( D sss (muj , m Hl , m q2 ) + D sss (m u , m H2 ,m qi )) 

~ ^{9l U ) 2 Dsss{m u ,m iL2 ,m iL2 ) 

+ 2(gf 2 U ) 2 D sss (m u ,m bL2 ,m bL2 ) 

+ (& D + g U S2 D ?D sss (m u1 m bm ,m- b J + 2{g u Si D f D sss {m^ m bm ,m bR2 ) 
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+ 2(gV D ) 2 D sss (m u ,m~ bR2 ,m~ b J] 
a 2 

- ^[2N C (N C - 2)D sv {m iL2 , 0) + (N c - l)[2D sv {m iL2 , m G ) + 2(c 2 q D sv (m< jl ,m G ) 

o 

+ s 2 q Dsv(m H2 ,m G ))} + ^[2D sv (m iL2 , m G ) + 2(c 2 q D sv (mq 1 ,m G ) + s 2 q D sv (m H2 , m G ))\ 

+ 2N C {N C - 2)D sv (m~ bL2 ,0) + (N c - l)[2D sv (m bL2 ,m G ) + 2(<? q D sv (m h ,m G ) 

+ s 2 q Dsv(m Hl ,m G ))} + ^[2D 5y (m^ 2 , m G ) + 2(c 2 D S v(m q2 ,m G ) + s 2 q D S v(m Hl ,m G ))] 

+ 2N C {N C - 2)D sv (m bR2 ,0) + (N c - l)[2D sv (m bR2 ,m G ) + 2D sv (m bm ,m G )} 

+ j^[ 2D sv{m bR2 ,m G ) +2Dsv(m bm ,m G )} 

~ ^9l 3 [4 D sv(m qi ,0) + D sv (m iL2 , 0) + %K H ,0) 

+ c 2 q Dsv(mq 2 , 0) + D S v(m bL2 , 0) + D S v(m bL3 , 0) + s 2 q D S v(m Hl , 0) + s 2 q D S v(m H2 , 0)] 

+ ^Ai + Aq 3 ^ [2cjD S s(m 5l , m^) + 2s q D ss (m h2 ,m h2 ) + 6s 2 q c 2 q D ss (m h2 ,m qi ) 

+ A(s 2 q D ss (m h2 ,m~ tL2 ) + c 2 q D ss (m qi ,m iL2 ) + 6D 5S (m t - i2 , m fi2 ) 

+ 2c q D ss (m q2 ,m q2 ) + 2s q D ss (m hl ,m hl ) + Qs 2 q c 2 q D ss (m hl ,m q2 ) 

+ A(s 2 q D ss (m hl ,m~ tL2 ) + c 2 q D ss {mq 2 ,m~ tL2 ) + 6£>ss(m gl . 2 ,m gl . 2 )] 

+ A D3 [8D ss (m 6iji , m bm ) + 2AD ss (m bR2 ,m bR2 ) + 16D ss (m gm , m bR2 )\ 

+ -j{h? U + 5'? 1 c/ )[cjL'55(m^ 1 ,m u ) + c 2 q D ss (m qi , mj + s 2 q D S s(m h2 ,m u ) + s 2 q D S s{m h2 ,mS) 

+ 2D55 (m^, fflj + 2D S s(m iL3 , mj)] 

+ ^f[cS%(mt., "la) + 5CgD55(m w , m^) + s 2 q D ss (m u , m h2 ) 

+ 5sjDss(m a; , m ft2 ) + D ss (m u , m~ tL2 ) 

+ 5L>55(m aJ , 77if M ) + D S s(m u , m t ~ i3 ) + 5Dss(mu, m t ~ i3 )] 

+ s 2 q D S s(m hl ,m u ) + s 2 q D S s(m hl ,m w ) + 2D ss {m L2 ^ m ^) + 2Dss(m t - J . 3 ,7n w )] 

+ ^^[Cg-Dssl^, m^ 2 ) + 5cjD ssK, m g - 2 ) + s 2 q D ss (m u , m hl ) 

+ 55^55(7710,, m hl ) + D ss {m u , m bL2 ) 
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+ ^Dssirriu, m bh2 ) + D ss (m u , m^J + 5.055(777^, m 6j . s )] 
+ + + 2^s5[D ss (m hl ,m ftl ) 

+ D S s(m h2 ,m h2 ) + D S s(m hl ,m h2 ) + Dssim^m^) + £55(777^, m^) 

+ D ss (mq 2 , m q2 )} + c 2 q [D ss (m iL2 ,m h2 ) + D ss (m iL3 ,m h2 ) 

+ D ss (m~ bL2 ,m hl ) + £55(777^,777^)] + sJ[D ss (7n fi . a , rn^-J + D 55 (m t - M , m qi ) 

+ D S s(m~ bL2 ,mq 2 ) + D ss (m bL3 ,m q2 )} + c 2 2q [D S s{m qi ,m h2 ) + D ss (mq 2 , m hl )\ 

+ ((A 3 + A^)(cJ + sj) - 2sJcJ(A^ + ^j^slm^mJ+^K.mJ] 

+ +2A^Cg[D5'5'(m^ 1 ,m ?l2 ) +^55(771^,771^)] 

- (a 3 + A^ [cj (D ss (m iM ,777^) + D ss (m~ bL3 , m fc2 ) 

+ D ss (rn- tL2 ,m hl ) + £55(777^,777^)) + s^r^m^m^) + £>5s(™g M , m^) 

+ D ss {m iL2 ,mq 2 ) + D ss {m iL3 , m q2 ))\ 

+ 2Ai(2 - JVcj^s^,-, m g - 2 ) 

+ s q D S s(m hl ,m h2 ) + s 2 q c 2 q D S s(m h2 ,m q2 ) + s 2 q c q D S s(m hl ,m qi ) 

+ D ss {m iL2 ,m~ bL2 ) + D ss (m iL3 ,m b J} 

+ 9^[c q D S s{m qi ,m q2 ) + s q D ss (m hl ,m h2 ) + sjcj (£^(771^,771^) 

+ D55 (m 9 - 2 , 777 h2 )) + D55 (m fi2 , 777^) + D55 (m^ , 777^)] 

+ ^[^ss^nWiffl) + s q D S s(m hl ,m h2 ) + s 2 q c 2 q (D S s(m qi ,m hl ) 

+ D ss (mq 2 ,m h2 )) + c 2 q D ss (m qi ,m bL2 ) + s 2 q D S s(m h2 ,m bL2 ) 

+ c 2 q D S s{m qi ,m bL3 ) + s 2 q D S s{m h2 ,m bL3 ) 

+ c 2 q D ss (m iL2 ,mq 2 ) + sjDss (77i fj . a , m hl ) + D ss (m iL2 ,m bL2 ) + D ss (m iL2 ,m bL3 ) 

+ c 2 q Dss{m~ tL;i ,mq 2 ) + s^L^m^, m hl ) + ^ss("if M , "ig M ) 

+ L>55 K i3 ,777- bM ) 

+ ^ [c^55 (m^, 777 Sm ) + s 2 q D ss (m h2 , m bm ) + D ss (m iL2 ,m bR2 ) 

+ D55 (m fi3 , 777 Sfl3 ) + cjD ss (m fe , m Sm ) + s 2 q D ss (m hl , m bm ) 

+ D ss (m bL2 ,m bR2 ) + D ss (m bL3 ,m b J] 
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+ 9% D {^ q D ss{mq X , m~ hm ) + s 2 q D ss (m h2 , m Sj J 

+ (? q D S s{m^m liu ) + s 2 q D ss {m h2 ,m lm ) + c 2 q D ss {m^m bR3 ) + s 2 q D S s{m h2 ,m- hR ,) 

+ D ss (m h2 , m- bRi ) + D ss (m h2 , m~ bR2 ) + D ss {m~ th2 , m~ bR3 ) 

+ D ss (m~ tL3 , m lm ) + D ss {m hz , m lm ) 

+ D 55 .(m fM ,m Sj{3 ) 

+ cjDssK, m 5j J + sJ-^Mm^, m 5j J 

+ cjD S5 (m f2 ,m Sfl2 ) +sjD ss (m hl ,m gfl2 ) + cjD ss (m fe ,m gji3 ) + sjD 55 (m ftl , m Bfl3 ) 

+ D ss (m~ bL2 ,m~ bm ) + D SS {m ha ,m- b J + D ss (m bL2 ,m b J + D ss (m bL3 ,m bm ) 

+ D ss (m bL3 ,m bR2 ) + D ss (m h3 ,m b J} 

+ g^^ 5 ^' m bm"> + D ss{m UJ , m bm ) + 2D SS (m u , m~ bRa ) + 2D S5 (m w , m^J] 

+ ^^[Are^u, m 6m ) + Dss(m u , m bR2 ) + D ss (m u , m bR3 ) 

+ 5D ss (m w , mi ) + 5£>s S (wu„ mt ) + hD ss {m u> , m b )]. 
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C Zero-temperature renormalization 

The most important zero-temperature renormalization effects with respect to our calculation con- 
cern the mass parameters. We will not go into the details of the renormalization, but refer the 
reader to the literature in which the pole masses for the relevant particles of our calculation have 
been obtained considering the full particle spectrum of the MSSM |37|, |38| . We use the expressions 
given in ref. |j37H , keeping only the top Yukawa coupling, in the appropriate (large- m^) limit. In 
this limit the relevant expression for the one-loop corrected Higgs mass is given by [|3^, |40|: 



m. 



167r 2 m 2 



2 

[(2 In h <2 + ' t2 



mi cos z 2(3 + 

' ' — ' II) J 

.2 \ 2\ -2 



nix 



4m? 



sin 2 2a t ) 2 f(m~ 1 ,m~ 2 ) 



2« 



m 



{)/(2m 2 t ) sin 2 2a t ln[m t - 2 /m t -J)] 



(71) 



where f(x,y) = 2~^\n^. 

In principle, there exists a metastable region in which the colour-breaking minimum is lower 
than the physical one at zero temperature. The constraint for absolute stability can be obtained 
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by studying the effective potential at zero temperature |2U|. This gives the constraint — mfj < 



(m^) 2 , where 



m u=[^f JL ) ■ (72) 
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Figure 1: Feynman diagrams contributing to the two- and four-point Green functions. 
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Figure 2: Critical temperatures in the 0- (solid) and x~ (dotted) directions as functions of m ; - 2 for 
tan j3 = 5 and rriQ = 300 GeV. 
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Figure 3: Plot of ^ as a function of in the 0- direction for tan (3 = 5, tuq = 300 GeV and 
A t = 100, 200 GeV. 
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103.0 




Figure 4: Allowed region in m^-m^ plane for vtiq = 300 GeV and = 200 GeV. To the left of the 
solid line there is a sufficiently strong first-order phase transition, to the right of the dotted line the 
physical vacuum is absolutely stable. The dashed line separates the region for which a two-stage 
phase transition can occur. 
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110.0 




Figure 5: Allowed region of parameter space in the m^-mf 2 plane for ijiq = 300 GeV, varying A t 
and tan (3. The dashed line is defined when the critical temperatures in the 0- and x~ directions are 
equal for the same variations of tan/3 and A t . 
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Figure 7: Allowed region of parameter space in the m^-m^ plane for itlq = 1 TeV, < A t < 650 
GeV and tan/3 = 5. To the left of the solid line there is a sufficiently strong first-order phase 
transition, to the right of the dotted line the physical vacuum is absolutely stable. The dashed line 
separates the region for which a two-stage phase transition can occur. 
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